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Abstract 

We study static and dynamical properties of fiuxons in a long annular Josephson junction (JJ) 
with a current injected at one point and collected back at a close point. Uniformly distributed 
dc bias current with density 7 is applied too. We demonstrate that, in the limit of the infinitely 
small size of the current dipole, the critical value of 7, above which static phase distributions 
do not exist, that was recently found (in the Praunhofer's analytical form) for the annular JJ 
with the length much smaller than the Josephson penetration length, is valid irrespective of the 
junction's length, including infinitely long JJs. In a long annular JJ, the dipole generates free 
fluxon(s) if 7 exceeds the critical value. For long JJs, we also find another critical value (in 
an analytical form too), which is always slightly smaller than the Praunhofer value, except for 
points where the dipole strength is iirN with integer N, and both values vanish. The static phase 
configuration which yields the new critical value is based on an unstable fluxon-antifluxon bound 
state, therefore it will probably not manifest itself in the usual (classical) regime. However, it 
provides for a dominating instanton configuration for tunnel birth of a free fluxon, hence it is 
expected to determine a quantum-birth threshold for fiuxons at ultra-low temperatures. We also 
consider the interaction of a free fluxon with the complex consisting of the current dipole and 
antifluxon pinned by it. A condition for suppression of the net interaction force, which makes the 
long JJ nearly homogeneous for the free fluxon, is obtained in an analytical form. The analytical 
results are compared with numerical simulations. The analysis presented in the paper is relevant 
to the recently proposed new experimental technique of inserting fluxons into annular Josephson 
junctions. 

PACS numbers: 74.50.+r, 85.25.-j, 03.65.Xp 
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INTRODUCTION 



Long Josephson junctions ( JJs) are well known to be a unique physical system that allows 
one to experimentally study dynamics of topological solitons in the form of fiuxons (alias 
Josephson vortices, each carrying a magnetic-flux quantum $o) 

. For both experimental 

and theoretical studies, the most convenient object is an annular (circular) long JJ, in 
which the net number of initially trapped fiuxons is conserved, hence new solitons may 
only be created as fiuxon-antifiuxon pairs ^, ^. Fluxon dynamics in the annular junctions 
manifests itself in the clearest way, as it is not complicated by reflections from boundaries. 
The interest to annular long JJs stems from the great potential they offer for fundamental 
studies of dynamical properties of solitons, such as, e.g., emission of the Cherenkov radiation 
BQ. A challenging p.oblen, fa the theory . ,nantnn, =oliton dynamo, m long JJ. winch 
has been recently observed in experiment at ultra- low temperatures . On the other hand, 
annular JJs offer applications in cryoelectronics, such as sources of highly coherent microwave 
radiation Q| and radiation detectors ^ . Besides that, long annular JJs have a potential for 
designing fluxon qubits 0, [l^ and fluxon ratchets IlJ 

A key problem in experiments with annular JJs is the trapping of a desired number 
of fiuxons. The only previously known method which made it possible to do that in a 
controllable way is rather complicated, relying upon the use of a scanning electron (or 
laser) microscope Other studies relied on fluxon trapping in the course of coo ling the 



irse oi cooling tne 



system below the critical temperature Tc of the superconducting electrodes 
However, the latter technique does not provide for reproducible results, and requires heating 
of the junction to high temperatures. Moreover, trapped fiuxons generated in such a manner 
get often pinned by Abrikosov vortices, which may be trapped in the electrodes in the course 
of cooling below T^. 

In a recent work a new method to insert fiuxons into annular JJs was proposed 
and demonstrated experimentally and numerically. It is based on injecting a relatively 
large current / into the junction locally; the current flows from the injection point into the 
superconducting electrode and also across the Josephson barrier, and is collected back at 
another point of the same electrode, which is separated from the injection spot by a small 
distance D. The schematic view of such J J is given in Fig. ^ In the following, we use 
the normalized notation for the distance between the current injectors, d = D/Xj, and the 
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current, e = I/Xj, where A j is the Josephson penetration depth. 

In this setting, one is actually dealing with a current dipole, the total current traversing 
the Josephson barrier being zero. The dipole gives rise to a local magnetic flux $ in the 
region between the injection and collection points. If the magnitude of this pinned flux 
attains the flux quantum $05 it may become energetically favorable (in the presence of an 
additional dc bias current uniformly distributed along the junction) to have the pinned flux 
compensated by a negative flux — $o- In fact, this implies that a fiuxon-antifiuxon pair is 
created in the annular JJ, so that the antifiuxon carrying the flux — $0 is pinned by the 
current dipole, while elsewhere in the long junction there appears a free fiuxon carrying 
the compensating flux $o- Further, if $ exceeds N^q, then N free fiuxons are expected to 
appear in the annular J J. 




FIG. 1: A sketch of the annular Josephson junction with a current dipole formed by local injection 
of current / into one of the superconducting electrodes. 



One objective of this work is to theoretically analyze the above problem and find a critical 
condition for the generation of free fiuxons by the current dipole in a long JJ biased by the 
uniformly distributed dc bias current 7. Results will be obtained in an exact analytical 
form (in section II) for the case when the dipole strength k, = ed is finite, while d ^ and 
e ^ 00. Two different solutions will be obtained. One of them yields a critical value of 7 
which coincides with the Fraunhofer expression (see below), that was recently found as a 
condition for the transition into a resistive state in a short annular JJ, whose full length was 
assumed to be much smaller than Aj jjj]. In our solution, the length of the ring junction may 
be arbitrarily large. Thus, we demonstrate that the Fraunhofer formula for the critical bias 
current is a universal one, which is valid irrespective of the dipole strength and junction's 
length. 
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The second analytical solution reported in this paper is exact only for long junctions. It 
yields the critical value of 7 which is slightly smaller than the one given by the first solution. 
However, the configuration on which the second solution is based (essentially, it is a fluxon- 
antifiuxon bound state) is, plausibly, always unstable. For this reason, the second solution 



may be irrelevant to the experimental technique outlined above in the usual 



classical) 
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regime of operation. However, precisely this second solution defines an instanton 
controls the tunneling rate for quantum birth of a free fluxon, in the same long JJ, at ultra- 
low temperatures. The instanton predicts that the exponentially small factor, which limits 
the quantum-birth rate, will scale as a square root of that for the long uniform junction, i.e., 
the quantum yield is strongly enhanced by the current dipole. Therefore, the second solution 
is physically meaningful too. The analytical predictions will be compared with results of 
direct numerical simulations of the critical condition (in the classical situation). 

A free fluxon moving in the annular JJ will periodically collide with the antifiuxon pinned 
by the current dipole and the dipole itself. An important issue is to minimize interaction 
between the free fluxon and the dipole- antifiuxon complex, so that the junction would seem 
as uniform as possible for the moving soliton. In section HI, we demonstrate, in an approx- 
imate analytical form, that, keeping the dipole size d and current e finite (or proceeding to 
the limit of the point- wise dipole), one may indeed select their values so that the interaction 
is strongly suppressed. This prediction, which has immediate relevance to the experiment, 
is also compared with direct numerical simulations. 

CRITICAL CONDITIONS FOR GENERATION OF FREE CLASSICAL AND 
QUANTUM FLUXONS 

General considerations 

The system outlined above is described by the following perturbed sine-Gordon (sG) equa- 
tion for the superconducting phase difference 0, which is written in the standard notation 

m 

(j)tt - (j)xx + sm(j) = 

-a<pt - 7 - £ [5 (x - d/2) ~5{x + rf/2)] , (1) 
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where S is the Dirac delta-function [it is used in Eq. (jT)) due to the assumption, which 
corresponds to the experimental situation, that the injection and collection leads have their 
width much smaller than Aj], 7 is the density of the uniformly distributed dc bias current, 
and a is a dissipative constant. In the limit of the point-like dipole {d 0), the limiting 
form of Eq. (Q) is 

(pit -4>xx + sin (j) = -a(j)t- 'J + k6'{x), (2) 

where k = ed is the dipole strength defined above, and 6' is the derivative of the delta- 
function. 

The model (|21) has been first studied almost two decades ago by Aslamazov and Gurovich 

n 

|19| . These authors considered interaction of fiuxons with an Abrikosov vortex trapped in 
one of the junction's electrodes, so that the vortex' normal core is parallel to the tunnel 
barrier. Later, the influence of the Abrikosov vortex on the fiuxon bound states localized 
near the junction boundary was studied by Fistul and Giuliani [2^. 

In the special case k = vr, Eq. (0) is also used as a model which describes composite 
long JJs including segments of 0-type and vr-type Josephson barriers [2^. Such composite 
junctions can be fabricated, for example, by using zigzag barriers between s-wave and d- 



wave superconductors [2^. This model gives rise to stable semi-fiuxons (vr-fiuxons with a 
semi-integer topological charge). As it was pointed out by Goldobin et al. j21|, semi-fluxons 
can also be created in a conventional long junction with local current injection in the form 
of a current dipole. 

The critical condition for the transition to a free moving fiuxon is defined in the following 
way: it is necessary to find a maximum value 7max of the bias current density 7 such that, 
for given k, static (time-independent) solutions to Eq. (j2I) exist for I7I < 7max, and do not 
exist if I7I exceeds 7max- As it was clearly shown by dint of direct simulations of Eq. (Q) 
reported in Ref. Q, the disappearance of the static solution means the appearance of freely 
moving fluxon(s) in the long JJ. 

Double integration of the static version {(pt = 0) of Eq. Q in a vicinity of the point a; = 
yields the following boundary conditions (b.c.) to be satisfied at this point: 

0(x = +0) - 0(x = -0) = k; = +0) = = -0) . (3) 

Off the point x = 0, the static solution obeys the equation 

-0" + sin0 + 7 = 0. (4) 



Further analysis of the static problem based on Eqs. Q and will be presented in the 
next subsection. 

In the sG system operating in the quantum regime (i.e., in a long J J kept at extremely 
low temperatures), fiuxon-antifiuxon pairs can be produced as a result of under-barrier 
tunneling j2^. The corresponding tunneling rate is determined by an instanton solution, 
which starts with some static field configuration (which is just a flat phase distribution, in 
the case of a homogenous system) and, going in imaginary time under the energy barrier, 
ends up on the physical shell, with a state consisting of far separated fiuxon and antifiuxon. 
In was shown that, in the presence of local inhomogeneities, the quantum birth of fiuxon- 
antifiuxon pairs may be facilitated (the corresponding tunneling rate being enhanced by an 
exponentially large factor) if one or both solitons appear on the mass shell (appear as real 
quasi-particles after completion of the tunneling) in a state pinned by an inhomogeneity 



2J|. In this connection, it is important to find a threshold 7thr (a minimum value of 7) 



past which there appears a state with a pinned fiuxon and free antifiuxon, i.e., an effective 
threshold for the quantum birth of fiuxons in the quantum regime. A difference from the 
problem of finding the above-mentioned value 7max, which determines the threshold for the 
creation of a free fiuxon in the classical regime, is that, on the way from 7 = to 7 = 
7max, the corresponding static solutions to Eq. (j21) need not be stable; just on the contrary, 
instanton solutions usually go through unstable states, such as fiuxon-antifiuxon pairs, as 
the evolution in imaginary time does not require stability in real time [23! [2^. 



Critical current in the annular junction of arbitrary length, the classical regime 

In the case of a finite-length annular JJ, the above-mentioned static problem takes the 
following form: one should find a solution to Eq. such that it satisfies b.c. which is 
tantamount to Eq. (jHl): 

(t){x = L)- <p{x = 0) = (p\x = L) = <p\x = 0) , (5) 

where L is the full length of the annular junction. Irrespective of the length, Eq. Q is 
equivalent to the Newton's equation of motion in "time" x for a particle with the coordinate 
0, in the presence of the potential 

[/(0) = -70 + COS0, (6) 
7 



hence the motion conserves the Hamiltonian, 

In the annular junction of any length, the values of H on both sides of the matching point, 
X = 0, must be equal, as they belong to one and the same solution. Besides that, due to 
the continuity of d(j)/dx at the matching point x = 0, see Eqs. 0, the values of {d(j)/dx)^ 
are also equal on both sides. Thus, equating the values of the Hamiltonian, we obtain a 
condition 

cos (01 + k) + 7 (01 + k) = cos 01 + 701 , (8) 

where 0i is the value on on one left side of the matching point, the value on the right side 
being 02 = 0i + /t, according to Eqs. 
A more convenient form of Eq. (jSJ is 

7 = [(1 — cos k) cos 01 + (sin k) sin 0i] . (9) 

One can now look for the largest possible value of the expression on the right-hand side of 
Eq. Q that can be obtained by varying 0i. This yields an expression of the Fraunhofer's 
type for the critical bias-current density, 

7c = 2K"Sin(K/2) , (10) 

which is attained at0i = 7r/2 — k/2. 

Thus, Eq. ()10j) gives the largest value of 7 beyond which no static solution may exist 
in a vicinity of the point dipole in the annular JJ of any length, implying a transition to a 
dynamical regime, i.e., generation of free fluxon(s) in the case of the long J J. Note, however, 
that it may happen that more than one critical values exist, then Eq. gives, according 
to its derivation from Eq. ©, only the largest one among them. In fact, it will be shown 
in the next subsection that (in the case of the long JJ) there indeed exists exactly one more 
critical value, which is smaller than (fTIH) . see Eq. (fT^ below. 

Quite naturally, in the limit k — Eq. (jlUj) yields 7^ = 1, which is the commonly known 
critical value of the bias current density in the long homogeneous J J Q.The fact that 7c 
vanishes at k = 27r is also easy to understand: the static equation (j3]) has obvious stable 
uniform solutions, 

0j,"^ = - sin-i 7 + 27rn (11) 
8 



(provided that I7I < 1), with an arbitrary integer n, hence the b.c. Q with k = 2tt imphes 
that two such solutions, (f)"^^ and ^g^'*, match to each other across the point a; = 0; in the 
annular system, they must be connected by a static 27r-kink, which is possible exactly at 
7 = 0. Similarly, in the case k = 2'kN with any integer N , Eq. (jlUj) again yields 7c = 0, 
which means that fluxons will appear spontaneously in this case, without the application 
the bias current. 

The result fllO|) does not depend on the length of the annular junction; for this reason, it 
coincides with the threshold for transition to a resistive state in a short annular JJ, which 
„a. tou„d . .ce. wo.. Q. « . Meva,. .o .en«o„ that the actual calculation of 
the threshold in that work was performed in an altogether different way, so that the sG 
equation was not used at all. 

The fluxon-birth threshold in the long annular junction, the quantum regime 

Coming back to the case of the long JJ, we notes that a solution for 0(x) must assume the 
background value (fTTj) (for instance, with n = 0) far from the dipole (formally, at x ±00). 
Thus, in the case of a very long (formally, infinite) junction, we need to find two solutions 
of Eq. (jlj), in the regions, respectively, x > and x < 0, so that one solution assumes the 
asymptotic value at x — > +cxd, and the other one assumes the same value at x — > —00. 
The solutions must be matched at the point x = 0, pursuant to b.c. Q, and it is necessary to 
find a maximum value of I7I for which the matching is possible with a given dipole strength 

K. 

Consideration of the equivalent mechanical problem mentioned above, i.e., motion in the 
potential ©, demonstrates that there are exactly two solutions to this problem. One of 
them yields the above result (fTn|) . A description of both solutions (for the case of a very 
long JJ) is given below. 

Both solutions employ the single function 0(x) which satisfies Eq. (0} and assumes the 
same asymptotic value (fTTj) simultaneously at x — * +00 and x —00. This function 
represents a fluxon-antifiuxon (ff) bound state; it is an even function, with a maximum 
value 0max at the central point. Using the conservation of the Hamiltonian ((Tj), one can find 
a relation between 0max and which takes the form 




(12) 
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In addition to the f f solution, Eq. (@)) also has a semi-divergent one, which starts with the 
value (fTTjl . say, at x = — oo, and diverges as ~ (7/2) at x ^ +00. This solution may 
also be employed to construct a full static state consisting of two pieces that are matched 
as per Eqs. Q at x = 0. It is easy to check that the state built as the combination of the ff 
and semi-divergent solutions disappears, with the increase of 7, exactly at the critical point 
(I10|l . so this is how the solution corresponding to the Fraunhofer's formula (|Tn|l looks in the 
long JJ. 

Besides this matched state, there is another one, specific to the long junction. It can 
be constructed from two pieces belonging to the ff solution [it is easy to see that no state 
satisfying the b.c. (jS)) can be constructed if both pieces belong to the semi-divergent solution; 
thus no extra solution is possible in addition to the two presently considered ones]. 

At the critical point corresponding to the disappearance of the new matched state, the 
b.c. (jHI) are satisfied in an extreme form: on one side of the point x = 0, we have the uniform 
background solution, = 0O) while on the other side a piece of the ff solution is to be used. 
Due to the continuity of the derivative 0'(x) at x = [see Eqs. Q] and the fact that 0' 
vanishes at the matching point in the critical case [as the derivative of 0(x) = 0o is zero], 
the central point (the one with = 0niax) of the ff solution, which is the only spot where the 
derivative vanishes, must be set exactly at x = when one attains the critical configuration, 
see Fig. El In view of the first b.c. (jHI), this means that the background value (fTT|) and 
0max are related (in the critical state) so that 0max = k + (j)Q = k — sin~^ 7. This expression 
should be inserted in Eq. ()12|) . After that, the "worst" transcendental terms 7sin~^7 in 
the ensuing equation mutually cancel, and the remaining equation can be solved in an exact 
form, to yield a new critical value of 7, which, this time, we call a threshold one: 



Y4sin^(fi;/2) + (k — sin k)^ 
As well as in the case of the Fraunhofer expression pO|) . the limit value of 7thr for k ^ 
is 1, which is the above-mentioned critical bias-current density in the uniform long J J. On 



solution takes the form of a bound state of a fiuxon and an antifiuxon separated by a large 
distance, which implies that even a small bias current is able to destroy this configuration 
pinned by the dipole. Accordingly, one finds from Eq. (|T^ that 



7thr = 



2sin2(K/2) 



(13) 




the other hand, if k is close to 27r, i.e., the phase jump in b.c. (jHl) is close to 27t, the ff 



7thr ~ (47r)-' {k - 27rf 



(14) 
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FIG. 2: A sketch of the critical field configuration that gives rise to the threshold value (|13j) of the 
bias-current density. 

when K, is close to 2tt, which should be compared to the asymptotic form of the Fraunhofer 
critical value in the same case: as it follows from Eq. (fTIH) . |7c| ^ \k — 27r| /8. On the contrary 
to this linearly vanishing expression, Eq. (fT^ shows that the newly found threshold value 
vanishes quadratically. As well as the Fraunhofer critical value, the threshold one ()14|) has 
higher-order minima at k = 2ttN with integer. 

The comparison of the two critical values (versus k) is shown below as a part of Fig. EJ^a). 
As is seen from the figure (and can be proved analytically), the threshold value ()14j) is always 
smaller (except for the points k = 271 N, where both critical values vanish), which implies 
that, in the application to the tunnel-birth problem for fiuxons in the quantum system, it 
will be a dominating one. In principle, for the same reason it could play a dominant role 
in the transition to the dynamical regime in the classical case, but, in fact, this is hindered 
by the instability of the ff static configuration which underlies Eq. ()13|1 . Unlike this, it 
was checked in direct simulations that the configuration based on the matched ff and semi- 
divergent solutions, which yields the Fraunhofer's expression (fTUll . is always stable in the 
long junction (for the case k = vr, it was also shown in Ref. |25i]). However, as it was 
explained above, the dynamical instability of the intermediate state does not bar using the 
configuration as a part of the under-barrier instanton trajectory. 

While detailed consideration of the under-barrier production of fiuxons in the quantum 
system is beyond the scope of this work, we briefly mention that the above-mentioned 
instanton may be constructed, following the pattern of the critical configuration that yields 
the threshold value (fT^ . as precisely a half of the fluxon-antifluxon bound state on one side of 
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the point x = 0, matched to the uniform state on its other side. This makes the integral 
of the system's action, taken along the under-barrier evolution in the imaginary time, which 
determines the exponential smallness of the tunneling rate equal to exactly half the 
integral in the absence of the current dipole. For this reason, the exponential smallness is 
expected to scale as the square root of that in the uniform system, so that the tunneling will 
be strongly enhanced. As for the direct role of the threshold value (fT^ . it implies that, if 7 
is taken larger than it, the tunneling mechanism will directly generate a free fiuxon in the 
system. 

SUPPRESSION OF THE INTERACTION BETWEEN A FREE FLUXON AND 
THE ANTIFLUXON PINNED BY THE CURRENT DIPOLE 

As it was mentioned in the Introduction, another problem relevant to the experiment is to 
suppress interaction between an already existing free fiuxon and its antifiuxon counterpart, 
which is pinned by the dipole pair of currents in the long JJ. A natural way to select 
parameters of the dipole configuration in Eq. (Q) with, generally speaking, finite d and e, so 
that to minimize the interaction, is to consider the case when the distance between the free 
fiuxon and pinned antifiuxon is large enough, so that the interaction between them may be 
treated by means of the perturbation theory |26|. In the framework of this approach, the 
free fiuxon with zero velocity, whose form is 

0fl = 4tan~^ [exp (cr(x — 0] > (15) 

where ^ is the coordinate of the fiuxon's center, and a = ±1 is its polarity, is regarded as a 
quasi-particle. Since the term in the Hamiltonian of the full sG model which corresponds 
to the current-dipole terms in the equation, is iJdip = ^ [<t>{.x = d/2) — (f){x = —d/2)], it is 
easy to obtain an effective potential of direct interaction of the free fiuxon with the current 
dipole, substituting the waveform (fT3j) in i^dip- Eventually, one can find an effective force of 
the direct fluxon-dipole interaction, 

Fdip = ^ = 2ae [sech - d/2) - sech + d/2)] 

^ 8ae sinh(rf/2) e"«, (16) 

where we made use of the assumption that the distance |^| of the fiuxon from the dipole is 
large. The polarity a of the free fiuxon and the sign of the current e mist correlate so that 
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the free fluxon is repelled by the dipole, while the corresponding antifluxon is attracted by 
the it (otherwise, the antifluxon cannot be in the pinned state). As it follows from Eq. (|TB|l. 
this implies that a = signe. 

Besides that, the fluxon also interacts with the antifluxon pinned by the dipole. A well- 
known perturbative expression for the interaction force between the fluxon and antifluxon 
is j3| -Pantifi ~ — 32e~^. Then, a simple condition for the effective suppression of the net 
interaction between the free fluxon and the complex including the current dipole and the 
antifluxon pinned by it may be formulated as a condition for mutual cancellation of the two 
forces, the repulsive one F^ip and the attraction force -Fantifl, which produces the following 
result: 



The meaning of this result is that the free fluxon is expected to move nearly as in the homo- 
geneous long J J if the parameters of the current dipole are selected according to Eq. (fTTjl . 
Note also that, in the limit of the point-like dipole, which was dealt with in the previous 
section, i.e., e oo and d ^ 0, the condition (|T7j) takes a very simple form, k = 8. This 
result is obtained in the framework of the approximations adopted above, but it is not very 
different from the exact result that follows from Eq. (jlUj) . i.e., n = 2tt (at this value of k, 
the fluxon can be very easily separated from the dipole- antifluxon complex). 

The validity of the prediction ()17|) was checked against direct numerical simulations of the 
full equation (^. For instance, in the case shown below in Fig. with e = 8 and d = 0.5, 
the product on the left-hand side of Eq. ()17j) takes the value 4.17, and the interaction with 
the dipole-antifluxon complex indeed gives rise to a very small perturbation in motion of 
the free fluxon. Simulation run for values of |£:| sinh {d/2) quite different from 4 (not shown 
here, as the pictures are rather messy) demonstrate a much stronger perturbation. 

NUMERICAL SIMULATIONS 

In order to verify analytical results for the fluxon injection, we performed numerical 
simulations by solving the full equation (^. In the simulations, each (5-function in Eq. (Q) 
was approximated by its smooth counterpart. 



e\ sinh (d/2) = 4. 



(17) 



e5{x) ~ ri9i^) = "/^sech^ C^^/O > 



(18) 
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such that Tj^ = e, which comphes with the definition S(x)dx = 1. This approximation 
imphes that the injected current is spread over the distance ~ ^Aj. We will present results 
obtained with the values £ = 1 and £ = 0.1 , which lie in a typical experimentally accessible 
range of this parameter [la] . The numerically calculated current- voltage characteristics will 
be shown in normalized units as 7(f), were v is the average fiuxon velocity normalized to the 
Swihart velocity c. With this normalization v = 1 corresponds to the asymptotic voltage of 
the single-fiuxon step. The simulations are performed with the normalized junction length 
L = 10 and dissipation coefficient a = 0.1. 

Figure Efa) presents the calculated dependence of the critical current 7c on the injection- 
current amplitude e for d = 2. The inset shows the adopted approximation for the injected- 
current profile, /(x) = ri[g {x — d/2) — g {x + d/2)], see Eq. (fT^ . with £ = 1 and t] = 9. As 
it has already been concluded from experimental data Q|, the dependence in Fig. |3]is very 
similar to the conventional Fraunhofer pattern of the critical current in a small Josephson 
junction in the magnetic field, the length of the equivalent small junction being associated 
with the distance d between the injecting points. The overlap between the lobes gets larger 
for larger d. Due to this overlap, the minimum value of the critical (fiuxon-depinning) 
current between the lobes decreases with d. 

Numerically calculated I-V curves for various values of the injection current (indicated in 
the plots) are shown in Fig. El^b). One- fiuxon and two-fiuxon steps can be clearly recognized 
here. Numerical data show that the steps on I-V curves are accounted for by fiuxons freely 
moving under the action of the uniformly distributed bias current. It can be noted from 
Fig. ISfa) that there is a residual pinning of fiuxon(s) due to the disturbance produced by 
the current injectors. This pinning is smallest at injection-current values which lie between 
the lobes of the 7c(£^) curve. Thus, for a given injector spacing d, the residual fiuxon pinning 
can be minimized by choosing an appropriate value for the injection current e. 

In order to numerically model the point- like dipole given by the last term of Eq. (j2)), we 
used the derivative g'{x) of the function g{x) defined by Eq. (fTH|l . Figure IH^a) shows the 
calculated dependence 7c on the dipole-current amplitude k. The inset shows the numerically 
used current profile, f{x) = rjg' {x — d/2) with £ = 0.1. The spatial grid used in this 
calculation had Ax = 0.01. In spite of the limited accuracy of our numerical scheme because 
of the sharp current profile, the agreement between the numerical data and the Fraunhofer 
formula (|TO|l . represented by the solid line, is very good. For the comparison's sake, in 
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FIG. 3: (a) Numerically calculated dependence of the critical current of the annular junction 
7c on the injection-current amplitude e for d = 2. (b) Numerically calculated current-voltage 
characteristics of the annular junction for two values of the injection current, e = 5.2 and e = 9.1. 

Fig. El^a) the dashed line shows the newly-found threshold value given by Eq. (jlHj) . 

The numerically calculated /-Incurves for the point-like dipole are shown in Fig. IH^b). 
In comparison with Fig. Efb), the depinning current is very small, and the value of the 
bias-current density corresponding to the steps is very close to the maximum value, 7 = 1. 
Both the one-fluxon and two-fluxon steps can be hardly distinguished from those in the ideal 
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uniform annular junction with trapped fluxons. 

The presentation of numerical results is concluded by Fig. It shows a 2D gray scale- 
plot of the spatially-temporal evolution of the instantaneous normalized voltage ipt{x,t) in 
the annular junction for e = 8, d = 1, and 7 = 0.4. The moving fluxon is recognized as 
a solitary- wave packet propagating at a nearly constant velocity across the junction. One 
can see that the disturbance of the fluxon motion in the region where the current injector is 
located is tiny, which was explained above, by the proximity of this case to the interaction- 
suppression condition (fT7|) . The residual perturbation may, nevertheless, be significant at 
small velocities of the fluxon, when its kinetic energy is comparable to the pinning potential 
generated by the perturbation. 



CONCLUSION 



In this work, we have considered static and dynamical solutions for fluxons in the model 
of the annular Josephson junction with the local current dipole. The analysis is relevant 
for interpreting the recent experiment that implemented a new technique of inserting 
fluxons into long junctions. 

Recently, a critical value of the bias-current density 7, above which the system per- 



forms transition to a dynamical state, was found |17|, in the Fraunhofer's form, for annular 
junctions whose length is much smaller than Aj. We have shown that, in the limit of the 
infinitely narrow current dipole, the same critical expression is valid at all values of the junc- 
tion's length. In the long junction, the dipole generates free fluxon(s) when 7 exceeds the 
critical value 7c. We have also found another critical value 7thr for long junctions, which is 
always slightly smaller than the Fraunhofer's one, except for points where both vanish. The 
phase configuration which yields the new critical value is generated by an unstable fluxon- 
antifluxon bound state, therefore it is not relevant to the classical-fluxon-generation regime. 
However, it determines a quantum-birth threshold for fluxons in the quantum regime. It was 
also concluded that these two critical values exhaust all possible solutions to the problem of 
finding the critical bias-current density in the presence of the current dipole. 

The interaction of a free fluxon with the complex consisting of the current dipole and 
antifluxon pinned by it was considered too. A condition for suppression of the effective 
interaction force was predicted, which makes the long junction effectively uniform for the 
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motion of the free fluxon. 

The analytical predictions obtained in this work were checked against direct simulations. 
In all the cases, they agree well. 
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FIG. 4: (a) Numerically calculated dependence of the critical current in the annular junction, 7c, 
on the dipole-current amplitude k. The newly found threshold value, given by Eq. (|13l (dashed 
line), and the Fraunhofer value, given by Eq. IjlUj) (solid line) are included too. (b) Numerically 
calculated current-voltage characteristics of the annular junction for two values of the dipole's 
strength, k = 6.3 and k = 12.8. 
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FIG. 5: Spatially-temporal evolution of the normalized instantaneous voltage in simulations of the 
annular Josephson junction with e = 8, d = 1 and 7 = 0.4. Dashed lines indicate the current 
injection points. 
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